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Abstract 
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CN . We study string theory duals of the expanding boost invariant conformal gauge the- 

O ■ 

CNj ■ ory plasmas at strong coupling. The dual supergravity background is constructed as 

an asymptotic late-time expansion, corresponding to equilibration of the gauge theory 
plasma. The absence of curvature singularities in the first few orders of the late-time 
expansion of the dual gravitational background unambiguously determines the equilib- 
rium equation of the state, and the shear viscosity of the gauge theory plasma. While 
the absence of the leading pole singularities in the gravitational curvature invariants 
at the third order in late-time expansion determines the relaxation time of the plasma, 
the subleading logarithmic singularity can not be canceled within a supergravity ap- 
proximation. Thus, a supergravity approximation to a dual description of the strongly 
coupled boost invariant expanding plasma is inconsistent. Nevertheless we find that 
the relaxation time determined from cancellation of pole singularities is quite robust. 
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1 Introduction 



Gauge/string correspondence [1, 2] developed into a valuable tool to study strong 
coupling dynamics of gauge theory plasma, and might be relevant in understanding 
the physics of quark-gluon plasma (QGP) produced at RHIC in heavy ion collisions 
[3, 4, 5, 6]. Besides applications devoted to determine the equation of state of the gauge 
theory plasma at strong coupling [7, 8, 9, 10, 11], the dual string theory models have 
been successful in computing plasma transport properties [12, 13, 14, 15, 16, 17, 18], 
the quenching of partonic jets [19, 20, 21], and photon and dilepton production [22]. 
Much less studied are the truly dynamical processes in plasma 1 . 

In [25] a framework of constructing the string theory duals to expanding boost 
invariant conformal plasma was proposed 2 . This was later used to study viscous prop- 
erties of the expanding plasma [26, 27, 28], other applications were further discussed in 
[30, 31, 32, 34, 33, 35, 36, 37]. The basic idea of the approach is to set up the asymp- 
totic geometry of the gravitational dual to that of the boost invariant frame, suggested 
by Bjorken [38] as a description of the central rapidity region in highly relativistic 
nucleus-nucleus collisions. One further has to specify the normalizable modes (a string 
theory duals to expectation values of appropriate gauge invariant operators in plasma) 
for the metric and other supergravity/string theory fields so that the singularity-free 
description of string theory is guaranteed. As explained in [25, 27, 28], requiring the 
absence of certain curvature singularities in dual gravitational backgrounds constrains 
the equilibrium equation of state of supersymmetric M = 4 Yang-Mills (SYM) plasma 
(in agreement with [7]), its shear viscosity (in agreement with [39]), and the relaxation 
time of M = 4 SYM plasma (in agreement with [40]). Despite an apparent success of 
the framework, it suffers a serious drawback: not all the singularities in the dual string 
theory description of the strongly coupled expanding M = 4 SYM plasma have been 
canceled [28] . The subject of this paper is to comment on the gravitational singularities 
observed in [28]. 

To begin with, we would like to emphasize the importance of singularities in string 
theory. As we already mentioned, constructing a string theory dual to expanding 
boost-invariant plasma implied identification of a set of string theory fields with non- 
vanishing normalizable modes, dual to vacuum expectation values (VEVs) of the gauge 

1 One notable example is the analysis of the shock waves generated by a heavy quark moving in 
gauge theory plasma [23, 24] . 

2 Earlicr work on more qualitative aspects of expanding plasma includes [29]. 
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invariant operators in plasma. A priori, it is difficult to determine which operators in 
plasma will develop a vacuum expectation value 3 . In practice, on the dual gravitational 
side of the correspondence one usually truncates the full string theory to low-energy 
type IIB supergravity, and often further to a consistent Kaluza-Klein truncation of 
type IIB supergravity, keeping supergravity modes invariant under the symmetries of 
the problem. From the gauge theory perspective, such an approximation implies that 
only certain operators are assumed to develop a VEV. Typically these operators are 
gauge-invariant operators of low dimension. Of course, the latter assumption might, 
or might not be correct: luckily, on the gravitational side the consistency of the ap- 
proximation is severely constrained by requiring an absence of singularities in the bulk 
of the gravitational background. A typical example of this phenomena is thermody- 
namics of cascading gauge theory [42]: if one assumes that the only operators that 
develop a VEV at high temperatures in this gauge theory are the same ones as those 
at zero temperature, all dual gravitational backgrounds will have a null singularity 
[43, 44]; on the other hand, turning on a VEV of an irrelevant dimension-6 operator 
[45] leads to a smooth geometry [46]. Thus, the presence of bulk curvature singularities 
in gravitational backgrounds of string theory points to an inconsistent truncation. 

In [25, 27, 28] the authors assumed a truncation of the full string theory dual 
to M = 4 expanding boost-invariant SYM plasma to a Kaluza-Klein reduction of 
type IIB supergravity on S* 5 , maintaining only the normalizable modes of the five- 
dimensional metric (dual to a stress-energy tensor of the plasma), and the dilaton 
normalizable mode (dual to a dimension-4 operator Tr F 2 ). The absence of singularities 
in five-dimensional Riemann tensor invar iant R% p R (5)fM/Xp of the asymptotic late- 
time expansion of the background geometry at leading and first two subleading orders 
correctly reproduced the equilibrium equation of state of M = 4 SYM plasma and its 
shear viscosity. Furthermore, removing the pole sing ularities in R ( ^ Xp R (5),luXp at third 
order in late-time expansion determined for the first time the relaxation time of the 
strongly coupled M = 4 SYM plasma [28] . Given the truncation of the string theory, 
it was not possible to cancel a logarithmic singularity at third order in R^ x ]ffl tlvXp . 
Turning on an appropriate dilaton mode removes a logarithmic singularity in the five- 
dimensional string frame metric in R^x"^ R( 5 > strm9 ^ uXp . However, as we show in the 
present paper, the full 10 dimensional metric remains singular even in the string frame. 

3 For a weak coupling analysis in QGP pointing to a development of a nonzero expectation value 
for TrF 2 see [41]. 
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The paper is organized as follows. In the next section we review the boost in- 
variant kinematics of gauge theory plasma. We identify parameters in the late-time 
expansion of the stress-energy tensor one-point correlation function with the hydro- 
dynamic parameters of Miiller-Israel-Stewart transient theory [47, 48]. Following the 
idea that a presence of singularities indicates an inconsistent truncation, we consider 
in the section 3 the full 5*0(6) invariant sector of type IIB supergravity, assuming 
the parity invariance of the M = 4 SYM expanding plasma. Compare to the analy- 
sis in [28], this includes one additional SO (6) invariant scalar, dual to a dimension-8 
operator. We show that using this additional scalar one can cancel the logarithmic 
singularity at the third order in late-time expansion either in ten-dimensional Einstein 
frame Ricci tensor squared TZ llu 7l flu , or in ten-dimensional Einstein frame Riemann 
tensor squared 1Z pu \ p TZ pl ' xp , but not both the singularities. Given the assumed symme- 
tries of the expanding plasma system, there are no additional supergravity fields to be 
turned on to cancel the singularity in string theory dual to expanding N = 4 SYM 
plasma. In section 4 we further show that canceling a singularity typically requires 
additional supergravity (or string theory) fields to be turned on. Using the consis- 
tent truncation presented in [45, 50], we consider supergravity dual to superconformal 
Klebanov-Witten gauge theory [51] plasma. Here, in addition to a dimension-8 opera- 
tor (analogous to the one discussed in section 3) there is an additional scalar, dual to a 
dimension-6 operator [45] . We show that turning on the latter scalar can remove loga- 
rithmic singularities in the third order in the late-time expansion of the ten-dimensional 
Einstein frame Ricci scalar 71, as well as IZ^IZ^ and lZ^ v \ p lZ pvXp . However, we find 
new logarithmic singularities at the third order in higher curvature invariants such as 
T^^iui\ipi'T^ tJ ' lUlX ' 2p2 T^ f j,2y2 XlPl 1Z p2V2 A 2 p 2 ) as wen as logarithmic singularities with differ- 
ent coefficients in (TZ....) 8 , (1Z. ...) 16 and so on. Canceling the logarithmic singularities 
in all the curvature invariants appears to require an infinite number of fields. Since 
such fields are absent in supergravity approximation to string theory dual of expanding 
boost invariant plasma, we conclude that such a supergravity approximation is incon- 
sistent. In section 5 we present further directions in analyzing the time-dependent 
framework proposed in [25, 27, 28]. 
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2 Review of boost-invariant kinematics 



Quark-gluon plasma is expected to be produced in high-energy collisions of heavy 
ions. In nucleus-nucleus experiments a central rapidity plateau structure for particle 
formation has been observed. Specifically, the expansion of the plasma appears to be 
longitudinal and homogeneous near the collision axis in the case of central collisions. 
Thus the system is expected to be boost-invariant in the the longitudinal plane. It is 
convenient to introduce the proper-time r and the rapidity y, defined as 

x° = Tcoshy, x 3 = Tsinhy. (2.1) 

Then rapidity invariance amounts to independence on y. We will also assume y — > — y 
invariance. 

We will also make another approximation, following [38], namely independence on 
the transverse coordinates x±. This corresponds to the limit of very large nuclei. 

With these assumptions, the whole energy-momentum tensor can be expressed 
purely in terms of a single function 4 e(r), the energy density in the local rest frame. 
This remaining function should be fixed by gauge theory dynamics. 

The proposal of [25] to use the AdS/CFT correspondence to determine e(r) amounts 
to constructing a dual geometry, by solving supergravity equations with the boundary 
conditions 

where the (5-dimensional part of the) metric is written in the Fefferman-Graham form: 

? 2 _ QfjLvdx^dx + dz , . 

as 5dim. — ^2 • ^ A - 6 ) 

Given any e(r) it is possible to construct such a metric. In [25] it was advocated 
that the requirement of the non-singularity of this geometry picks out a physically 
acceptable e(r). In [25, 27, 28] this problem was analyzed in an expansion for large 
t. At each order the square of the Riemann curvature tensor had poles which could 
be canceled only by a specific choice of the coefficients of e(r). In the third coefficient 
there remained a leftover logarithmic singularity in the Einstein frame 5 . In this paper 
we would like to revisit this issue. 



4 After taking into account energy-momentum conservation and tracelessness. For explicit formulas 
see [25]. 

5 The is also a singularity in the string frame Ricci tensor squared Ji^' strm 3) ft(w, string)^ _ 
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The outcome of this calculation is a specific r dependence for e(r): 

e(r) — (l — • • ■■) ■ (2.4) 

Now we may obtain the physical interpretation of the above result by checking 
whether it is a solution of some specific phenomenological equations. In the above 
case, the leading behavior corresponds to perfect fluid hydrodynamics, the subleading 
one is related to the effect of shear viscosity, while the third one is related to second 
order viscous hydrodynamics relaxation time. Let us note that the relaxation time 
depends on the form of phenomenological equations used to interpret (2.4). 

Muller-Israel- Stewart theory [47, 48] of dissipative processes provides a nice frame- 
work to study relativistic fluid dynamics. In the Bjorken regime of boost invariant 
expansion of plasma, the transport equations take the form [49] 

dr r r , . 

d$ $ 1/1 Id //3 2 \\ 211 



° dr + r n + 2 $ \t + (3 2 T dr \T J J 3 far' 

where T n is the relaxation time, e is the energy density, p is the pressure, $ is related 
to the dissipative part of the energy-momentum, and 

A = ^. < 2 -6) 

with r/ being the shear viscosity. In (2.5) anticipating application to conformal gauge 
theory plasma we set the bulk viscosity ( to zero. 
In the case of the M = 4 SYM plasma we have 

e ( r ) = -vr 2 iV 2 T(r) 4 , p(r) = -e(r) , v (t) = A s(r) = A -tt 2 N 2 T(r) 3 , 



T w (t) = r r Boltzmann( r 



Mr) 
Mr) 



(2.7) 



Given (2.7), we can solve (2.5) perturbatively as r — > oo: 



T(t) 



A [ i + h 1 $(r) = —ir 2 N 2 A — [ 1 + V /fc 

r i/s ^ + (Ar2 /3 )fcy | . 3- ^ r2 ^ + (AT 2/s) 



(2.8) 
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where A is an arbitrary scale. The constant coefficients {t^, fk} are related to the M = 4 
SYM dimensionless transport parameters A ( the ratio of shear-to-entropy density ) 
and r ( the plasma relaxation time in units of Boltzmann relaxation time ). For the 
first few coefficients we find: 

*i = -j^> t 2 = ~A 2 r, t 3 = _|-A 3 r(i + 24r), 

h = 2A(2r - 1) , f 2 = V(l - 8r + 24r 2 ) , (2.9) 

o 

f 3 = — A 3 (-l + 24r - 234r 2 + 1296r 3 ) . 

From the supergravity computations first done in [28] (and reproduced in the fol- 
lowing section of this paper) we find 

Matching the gauge theory expansion for the energy density (2.7) with that of the dual 
gravitational description (2.10) we find 



\/2 3 3 / 4 11 1 C 

A= 3V%' A = ¥JH T]0 > r = -18- 108 (2 - U) 

Using the supergravity results (3.29) and (3.33), we find 



A = ±-, r=i(l-ln2). (2.12) 

A different formulation of second order hydrodynamics in [52] derived from Boltz- 
mann equations lead to the equation 

0=^ + ^-1* 

_d$ $ 2 11 ( 2 - 13 ) 
°~^ + ^) ~ 3&r' 

which provides a different definition of the relaxation time ri n \ Using (2.13) leads 
to = |(1 - In 2) quoted in [28]. Of course (2.5) and (2.13) provide different 
interpretations of the same e(r). Without lifting the symmetry assumptions of the 
uniform boost invariant plasma expansion we may not rule out one of the (2.5) and 
(2.13) descriptions. 



7 



3 M = 4 QGP 



Supergravity dual to boost invariant expanding Af = 4 QGP was discussed in [25, 
27, 28]. We extend the previous analysis including all supergravity modes, invariant 
under the symmetries of the problem. Specifically, we assume that the SO (6) R- 
symmetry of M = 4 SYM is unbroken. We further assume that the boost invariant 
plasma is invariant under separate reflections of all three spatial directions. The latter 
assumption in particular ensures that the supergravity dual excludes the axion and 
various 3-form fluxes along extended spacial directions of the boundary — say RR 
fluxes with components F TX1X2 within the metric ansatz (3.11). In section 5 we comment 
why we believe that even allowing for a parity violating supergravity modes would not 
change our main conclusion, i.e., the inconsistency of the supergravity approximation 
in dual description of boost invariants expanding plasmas. 

Given the symmetry assumptions of the problem as well as the approximation of 
the full string theory by its low energy type IIB supergravity, the non-vanishing fields 
we need to keep is the five-dimensional metric, the warp factor of the five-sphere, the 
self-dual five-form flux and the dilaton. Parity invariance and the symmetries of the 
boost invariant frame further restricts the five- dimensional metric as discussed below. 

3.1 Consistent Kaluza-Klein Reduction 

Consider Einstein frame type IIB low-energy effective action in 10-dimensions 

Si ^iI/°^{ K -^) 2 -iV4- (31) 

with metric ansatz: 

dS?o = ~9MNdi M di N = 

= a- 2 (x) gflu (x) dx»dx u + <r 6/5 (x) (dS 5 ) 2 , 

where M, N, . . . = 0, . . . , 9 and fi, u, . . . = 0, . . . , 4, and (dS 5 ) 2 is the line element 
for a 5-dimensional sphere with unit radius. For the 5-form F 5 we assume that 

F 5 = JS + *.F 5 , -F 5 = -4Q tu s s , (3.3) 

where 0035 is the 5-sphere volume form and Q is a constant. We further assume that 
the dilaton is = (f>(x). 
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With the ansatz (3.2), we find 

U = a 2 ji? - ^ (d In a) 2 } + 20^ 

' F 2 = -8Q 2 a- 6 



(3.4) 



4-5! 



where i? is the Ricci scalar for the 5-dimensional metric g^. The 5-dimensional effective 
action therefore takes the following form: 

Sf = d'x ^g^R- l -{d4>f - ^{dlnaf -V{a)Y (3-5) 



"5 J Ms 

with 



«? 



P(<t) = -20a- 16 / 5 + 8QV- 8 , «g = . (3.6) 

It is convenient to introduce a scalar field a(x) defined as 

a(x) = e a{x) , (3.7) 

so that (3.5) can be rewritten as: 

Sf = d 5 x^g[R-\{d4>f - ^{daf -V{a)Y (3.8) 

where 

V(a) = -20 e- 16a/5 + 8Q 2 e~ 8a . (3.9) 

From (3.8), the Einstein equations and the equation of motion for a are respectively 
given by: 

24 1 1 

Rfiu = — (d^a) (d u a) + - (d^) (d v <f>) + -g^V(a) 



_5 dV 
4 

□0 = 



Ua = j-jf- (3-10) 
48 da 
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3.2 Equations of motion 

For the five-dimensional metric we use the same ansatz as in [25, 27, 28]: 
ds 2 = g^dx^dx" = 

J_ [- e Mr,z) dT 2 + e 2b(r,z) r 2 dy 2 + ^(r.z)^] + d 



(3.1V 



z* z* 

where dx\ = dx\ + dx\. It is easy to see that this is the most general boost invariant 
geometry written in Fefferman- Graham coordinates subject to a parity invariance along 
the boost direction 6 . 

Further assuming a = a(r,z) and = <f>(r,z), the equations of motion for the 
background (3.11), the dilaton and the scalar a become: 

• Einstein equations 



„2a 



d 2 z a + (d z a) 2 + (d z a) (d z b) + 2 (d z a) (d z c) d z a d z b d z c + 



d 2 b + {d T bf - (d T a) (d T b) + 2 (d 2 c - {d T cf - {d T a) (d T c)) - -d T a + -d T b 

T T 



,2 , 1 /Q i , 2 _le 2a 
3~z 2 



24 

= _ (d T ay + ^ (MY - ~V(a) , 



(3.12) 



,2a 



d 2 z b + (d z b) 2 + (d z a) (d z b) + 2 (9,6) (d z c) - -d z a - -d z b - -d z c + 



d 2 T b + (d T b) 2 - (d T a) (d T b) + 2 (d T b) (d T c) - -d T a + -d T b + -d T c 

T T T 



le 



2a 



3 z 



2 P(a) , 



(3.13) 



'2a 



d 2 c + 2 (d z c) 2 + (9 z a) (d z c) + (9,6) (0,c) - -d z a - -9,6 - -<9 z c + - 



d 2 c + 2 (d T cY - (d T a) (d T c) + (d T b) (9 T c) + - (9 T c) 

r 



(3.14) 



6 The parity invariance excludes metric components g yz (j^ z) and g T y{T, z), in principle allowed by 
a Fefferman- Graham coordinate frame. 
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d 2 z a + (d z a) 2 + d 2 z b + (d z b) 2 + 2 [d 2 z c + (d z c) 2 ] - -d z a - -d z b - -d z c + 1 

/C /C 



■y(^) --(9 2 0) 



9 T 9 2 6 - (9 T 6) (d z a) + (9 T 6) (9,6) + 2 [<9 r «9 2 c + (<9 r c) (9,c) - (d T c) (d z a)\ 

1 1 24 1 

d z a + -d z b = — - (d T a) {d z a) - - (d T <fi) (d g <f>) ; 

T T Z 

• dilaton equation: 

= z 2 ld 2 z <P + (d z a + d z b + 2d z c) (d z <f>) - - z (d z <f>) 



(3.15) 



(3.16) 



— e 



-2a 



<9 2 + ( -d T a + <9 r 6 + 2d T c + - ) (d T </>) 



(3.17) 



• equation of motion for a 



5 <9P 



48 fa = ^ * a + (<9z ° + dzh + 29zC) ^ za) ~ ~ ^ za) 



— e 



-2a 



i);<\ - ( i) : a ■ i) b ■ ■!<);<■ ■ - ) {0 r <\) 



(3.18) 



3.3 Late-time expansion 

In [25, 26, 27, 28] equations similar to (3.12)-(3.17) were solved asymptotically as a late- 
time expansion in powers of (t~ 2 ^ 3 ) (or powers of (t -1 ^ 3 ) for the dilaton) introducing 
the scaling variable 

v = 4^ (3.19) 



T 



1/3 



in the limit r — > oo, with t> kept fixed. Such a scaling is well motivated on physical 
grounds as in this case we find that the boundary energy density e(r) extracted from 
the one-point correction function of the boundary stress-energy tensor 



c(r) 



N 2 2a(2,T) 
hm 



N 2 , 2a(u,r) 
hm 



(3.20) 



271-2 "^5 ^4 2 7T 2 u-0 t;4 r 4/3 

would have a late-r expansion appropriate for a conformal plasma in Bjorken regime 
[28]. Thus we expect 

cl{t,v) = a (v) + -^a x {v) + ^a 2 {v) + -^a 3 {v) + C(r" 8 / 3 ) 



T 



T 



4/3 



6(t,t;) = 6d(«) + ^M") + ^M") + + 0(^ 8/3 ) (3.21) 

c(t,v) = cq(v) + ^aCi(u) + ^c 2 (v) + -^(u) + 0(r^ 3 ) . 
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In the equilibrium, both the dilaton and the a-scalar vanish, implying that in (3.12)- 
(3.16) they enter quadratically for small 0, a. Thus, given (3.21) we have two possible 
late-time asymptotic expansions for <f> and a: 

or (3.22) 

<Kt,v) = ^j-M v ) + \fa v ) + ^rM v ) + o{t~ 7 / z ) , 



and 



a ( r > v ) = -^ a M + ^ a ^ v ) + + °( T 8/3 ) 



or 



(3.23) 



a(r,v) = ^diW + ^a 2 (v) + ^& z {v) + 0(r 7/3 ) . 



For each of the four possible asymptotic expansions, we solve (3.12)-(3.18) subject to 
the boundary conditions 

Oi(v),bi(v),Ci(v)> =0, l(pi{y) oi (f)i{y)\ =0, 

aj(f ) or dj(f ) > = . 

Additionally, we will require the absence of singularities in the asymptotic late-time 
expansion of the quadratic curvature invariant X^: 

i m = n^xW'^ 



^PW + ^W^r'" - - <X25) 

Of course, ultimately, all the curvature invariants of the bulk geometry must be nonsin- 
gular for the supergravity approximation dual to boost invariant M = 4 SYM plasma 
to be consistent, unless singularities are hidden behind the horizon. However, it turns 
out that asymptotic expansions (3.21)-(3.23) are fixed unambiguously, given (3.24) 
and (3.25). Thus non-singularity of the other curvature invariants provides a strong 
consistency check on the validity of the supergravity approximation. 

3.4 Solution of the late-time series and curvature singularities 

[21 

The absence of curvature singularities in Xq determines the leading solution to be [25] 
ao = ^ln ( 1 ^ 4 4 / / 3 3 )2 , 6 = co = iln(l + ,V3). (3.26) 
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Notice that v = 3 1 / 4 is the horizon to leading order in r. It is straightforward to find 
that the absence of naked singularities in the bulk (specifically at v — 3 1 / 4 ) in zj 2 2 
would require that a\ = a 2 = &\ = &2 = 0, <f>i = 02 and 0i = 2 = 0. We further 
have [26, 27, 28] 



(9 — v 4 ) v 4 v 4 i]o 3 — f 4 

9 — v 8 3 + v 4 2 3 + v 4 

v 4 



and 



_ (9 + 5t; 4 )t; 2 (9 + v 4 )v 4 2 (-1053 - 171^ 4 + 9t; 8 + 7v 12 )v 4 
° 2 ~ 12(9 -v 8 ) 72(9 - t; 8 ) + % 6(9 - 1; 8 ) 2 



1 , v^-^ 2 3 2 3 



4 



+ 



^2 



24v/3 V v^ + t^ 2 4(3 + t; 4 ) 2 ^ (v 7 ^^ 2 ) 2 

2 - 4 (39 + 7u 4 )u 4 1 . VS-v 2 



v „ f 



& 2 = - 2c 2 + 777-T7TV + CtTTTTT— IV + ^0 /o I ..4 2 + 7w^ ln 



4(3 + v 4 ) 24(3 + v 4 ) ,u 2(3 + v 4 ) 2 8^ + v 2 



(3.27) 



8^ V3 + v 2 4 /u 3 + v 4 

7T 2 v 2 (9 + v 4 ) v 4 2 (-9 + 54v 4 + 7v 8 )v 4 

° 2 ~ 2887! + 12(9 - v 8 ) + 72(3 + t; 4 ) " % 6(3 + t; 4 )(9 - v 8 ) 

1 , v^- v 2 1 nn 2 , , 3-w 4 

H f= In -7= 1 (C + 66r7 2 ) ln (3 28 

8^ VS + v 2 72 1 l0) 3 + v 4 



3 2l 3-w 4 

H — ?? n In - , 

4 /0 3 + v 4 ' 

with [27] 

»*> =21^374- 

While the equations at the third order for {03,63,03} are a bit complicated, it is 
possible to decouple the equation for a 3 . The resulting equation is too long to be 
presented here. It is straightforward to solve the equation perturbatively as 

y = 3-v 4 ^0 + , (3.30) 
da 3 VQ 4 /2V231/4 2 1 / 2 3 1 / 4 hi2\ 3 ^, , nns 

which is sufficient to determine the singularities in xf^ at v — 3 1//4 . 
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The absence of pole singularities in Z 3 2 ' at v = 3 1//4 implies that only o 3 7^ 0: 

as = «3o(Y J ^ + — ) In — 4 -—) , (3.32) 
' V V 96w 4 864 / 3 - v 4 144 / ' v ; 



where 0:3,0 is an arbitrary constant; it further constraints [28] 

17_ 

7! 



C = 2V3 ln2 - (3.33) 



Using (3.31)-(3.33) we find 

if = finite + 2 1 ' 2 3 3 / 4 + y a 3 ,o) h(3 - ^ 4 ) , u - 3 V4 . (3.34) 

Notice that with o 3i o = the logarithmic singularity agrees with the one found in [28]. 
From (3.34) it appears that the curvature singularities in the bulk of the supergravity 
dual to expanding J\f = 4 SYM plasma can be canceled for an appropriate choice of 
0:3,0- While the Ricci scalar 1Z is indeed nonsingular, unfortunately, this is not the case 
for the square of the Ricci tensor. We find 



n^-RT = finite + \ o 3 , hi(3 - v 4 ) , v -> 3 V4 . (3.35) 



Thus, it is impossible to cancel logarithmic singularity both in r R. lLV W LV and TZ^ V p\TZ^ upX 
at order O (r -2 ) with a supergravity field o. 

3.5 Curvature singularities of the string frame metric 

In [28] it was suggested that the logarithmic singularity in the curvature invariants of 
the Einstein frame metric (3.2) might be canceled in string frame metric 

g strin g = ^ g ^ (33^ 

Unfortunately, this is not possible. Indeed, as in [28], in order to avoid pole singular- 
ities in the curvature invariants of the string frame metric at leading and first three 
subleading orders, the dilaton can contribute only at order 0{t~ 2 ) 

0(r, v) = \ k 3 In + ° (^ 8/3 ) > ( 3 - 37 ) 

where fc 3 is an arbitrary constant. 
We find then that 

R (10,string) R (10,strin 9 )^ = finite + J_ ^ ^ _ m fc ^ ^(3 _ ^4) ? (ggg) 
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and 7 

R (10,strin 9 ) R{W ,strin g )^ P X = finite + 1 ^ 2^ 3 3/4 + y a 3 ,0 - 152 k^j ln(3 - V 4 ) , 

(3.39) 

as v — + 3i/ 4 . While it is possible to remove logarithmic singularities in (3.38) and (3.39) 
by properly adjusting 0:3,0 and /C3, the higher curvature invariants of the string frame 
metric would remain necessarily singular. Specifically, for string frame fourth order 
curvature invariants 2^ 10 > strm 3)[4] we find 

-y-(10, string) [4] _ p(10, string) r>{\0, string) uiui td(10, string) fiu p(10, string)ti2i'2p^ 

=finite + 1 2V2 3 3 / 4 + a 3 ,o - 3904 fc s ) ln(3 - v 4 ) , ^ 

as v — > 3_ . 



4 KW QGP 

In the previous section we studied supergravity dual to Af = 4 SYM plasma assum- 
ing unbroken 5*0(6) global symmetry and the parity invariance along the extended 
boundary spacial directions. At a technical level, we found that the absence of bulk 
singularities in the gravitational dual of expanding boost invariant plasma is linked to 
nontrivial profiles of massive supergravity modes (corresponding to VEV's of irrele- 
vant gauge invariant operators in plasma). This suggests that turning on additional 
massive supergravity modes might remove curvature singularities. Within the assumed 
unbroken symmetries of the J\f = 4 plasma there are no additional massive modes in 
the supergravity approximation. However, we can test the link between massive super- 
gravity modes and singularities in the supergravity dual of expanding boost invariant 
conformal plasmas for a slightly more complicated example. 

In this section we study supergravity dual to boost invariant expanding Klebanov- 
Witten super conformal plasma [51]. Assuming that parity along extended boundary 
spacial directions as well as the global SU (2) x SU (2)xU(l) symmetry of KW plasma at 
equilibrium is unbroken, the effective supergravity description will contain two massive 
modes [45, 50] — a supergravity mode dual to a dimension-8 operator ( an analog of a 
field in the context of M = 4 plasma ) and a supergravity mode dual to a dimension-6 
operator. We expect that appropriately exciting both modes we can remove curvature 

7 The result for the Riemann tensor squared here corrects the expression presented in [28]. 
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singularities in all quadratic invariants of the metric curvature at the third order in 
late-time expansion. While we show that the latter expectation is correct, we also find 
that higher curvature invariants will remain singular in this model. In fact, it appears 
we need an infinite set of massive fields (which is not possible in the supergravity 
approximation) to have a nonsingular metric. 

Since the computations for the most part mimic the analysis of the previous section, 
we highlight only the main results. 

4.1 Consistent Kaluza-Klein reduction 

Consistent KK reduction of the KW gauge theory plasma has been constructed in 



[45, 50]. 

The five dimensional effective action is [50] 




where we defined 



V = -24e _ 3 f - 



+ 4e 3 




(4.2) 



We set the asymptotic AdS radius to one, which corresponds to setting 



K = 4. 



(4.3) 



From Eq. (4.1) we obtain the following equations of motion 




(4.4) 




(4.6) 



(4.5) 




(4.7) 



The uplifted ten dimensional metric takes form 



2 



ds 2 10 = gMdyV + n?(y)ej + Q 2 2 (y) £ {e\ + eJJ , 



(4.8) 



a=l 
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where y denotes the coordinates of A4$ (Greek indices /i, v will run from to 4) and 
the one-forms e^, eg a , (a — 1, 2) are given by 

e -4> = \ {d^ + cos da d ^ ' ee - = ~^ d9a ' 60a = "71 9(1 d ^ a ' 
Furthermore, 

gMdjfdy" = n- 2/3 n 2 s/3 ds 2 , n t = e f~ 4w , n 2 = e f +w , (4.io) 

and ds 2 is the five-dimensional metric (3.11). 
4.2 Late-time expansion and solution 

We consider the same five-dimensional metric ansatz as in (3.11); we use the late-r 
expansion of the metric warp factors {a,b,c} as in (3.21). In order to avoid pole sin- 
gularities in curvature invariants, the dilaton must be set to zero. For the supergravity 
scalar f(r,v) dual to a dimension-8 operator in KW plasma and for the supergravity 
scalar w(t,v) dual to a dimension-6 operator in KW plasma we use the asymptotics 

f(r,v) = ^f 3 (v) + 0(r-^) 



(4.11) 



w(t,v) = ^w 3 (v) + O (r- 8 / 3 ) . 



Since the massive supergravity modes {f,w} are turned on only at order 0(r~ 2 ), the 
metric warp factors {aj(i>) , bi(v) , q(i>)} are exactly the same as for the M = 4 SYM 
plasma, see (3.26)-(3.28), (3.29), (3.31) and (3.33). Moreover, analogously to (3.32) 

/aW=/3o(Y J ^ + — ) In^±4- J-V (4.12) 

For w 3 we find the following equation 

„ 5v 8 + 27 ,12 , 
= w» 3 + 8 «; - - ^ 3 • 4.13 
v (v 8 — 9) f 1 

Solution to w 3 must have a vanishing non-normalizable mode as v — > 0. Near the 
horizon (3.30), the most general solution to (4.13) takes form 

^3 = ^0,0 + «x>,i ln ( 4 33/4 2/) + ° (y) • ( 4 - 14 ) 



17 



We find that vanishing of the non-normalizable mode of w 3 as v — > requires 



in which case, 



Wo,o = Wo,i x 



w 3 



7Y 



QV3 



x 



wo,! v 6 + (v 14 ) 



(4.15) 



(4.16) 



4.3 Quadratic curvature invariants of (4.8) 

We collect here results for 71, r R. iW W w and H. iwp \Ry ,vpx curvature invariants of the 
metric (4.8) to the third order in the late-r expansion. 



4-3.1 At leading order 



K 



(o) 



-20 + 20 , 



(0) 



n^ pX n pupX 



(0) 



= 80 + 80 , 
8(5w 16 + mv 12 + 1566w 8 + 540w 4 + 405) 



+ 136, 



(3 + W 4 ) 4 

where in (4.17)-(4.19) we separated the AdS$ and the T 1,1 contributions. 



4.3.2 At first order 



K 



(i) 



= 0, 

(1) 



= 0, 



4-3.3 At second order 



(i) 



1 41472(v 4 - 3> s 



r 2/3 (3 + v 4y 



K 



(2) 



= 0, 

(2) 



0, 



(4.17) 
(4.18) 
(4.19) 



(4.20) 
(4.21) 
(4.22) 



(4.23) 
(4.24) 
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(2) 



1 



576<y - 3y 



c 



r 2/3 ^ (3 + V 4 )5 

6912(5v 24 - 60v 20 + 2313v 16 - 6912v 12 + 26487v s 



(3 - t; 4 ) 4 (3 + v 4 ) 6 
4608(5w 16 + 6w 12 + 162w 8 + 54w 4 + 405)w 10 



+ 

4608(5w 16 + 

(3 - w 4 ) 4 (3 + v 4 ) 5 
4.3.4 At third order 

Provided that C is chosen as in (3.33), we find that 



18468t; 4 + 13851) A 2 (4 25 ) 



K 



(3) 



0, 



once we use Einstein equations. Also, it is easy to determine that 

(3) mn // n \ q i . ^4 

J 3 ~ 

The non-singularity condition therefore requires that 



100 

"27 



/30 = 0. 

Finally, using (3.31), (4.14) and (4.28), we find 

(3) 



Thus, choosing 



(8 2 1 / 2 3 3 / 4 - 384u> ,i) ln(3 - v 4 ) 

2 - 2 1 / 2 3 3 / 4 In (6 (3c) 1 / 5 ) + 384u> , ) + O ((3 - v 4 )) 



W ,i 



2 l/2 33/4 

48 

the logarithmic singularity in (4.29) is removed. 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



4.4 Higher order curvature invariants of (4.8) 

Let us define shorthand notation for the contractions of the Riemann tensor. For each 
integer n we have 

^ fiup\ = ^ piuifiu ° ^ ^ p\ 1 (^-31) 

where 

pupX = H^ p x . (4.32) 
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We further define higher curvature invariants X^ n \ generalizing (3.25): 

= iPiv) + -ijifc-) + ^zFw + ^xT'w + of-- 8 ' 3 ) ■ (4 ' 33) 

With a straightforward albeit tedious computation we can extract logarithmic singu- 
larities in xf ' . For the first couple invariants we find: 

- 384 (w 0>1 - 3 3/4 V^j In (3 - v 4 ) + finite 

- 3072 (w ,i - ^ 3 3/4 V^j In (3 - v A ) + finite 
if = - 98304 (w 0;1 - ^ 3 3/4 V^j In (3 - v A ) + finite 

- 50331648 (w ,i + ^ 3 3 / 4 ^ In (3 - v 4 ) + finite ( 434 ) 

(1375731712 \ 
Wo,i + 27 33/4 ^ J ln ( 3 ~ + finite 

^641 rr .^ Mnmniinrninrr , inrn ( 5044031582654955520 ~ /4 ^ 
X [ p = - 56668397794435742564352 I w ,i + 3 3/4 V2 



T [2] 
x 3 

x 3 



-[16] 
3 

-[32] 



9 

x ln (3 - v A ) + finite , 

as v — > 3^ 4 . Interestingly, we find that all lower order invariants, ie., 

if 1 , n = {1,2, 3, 4, 5, 6}, 2 = {0,1, 2}, (4.35) 

are finite as v — > 3 1 / 4 . 

Clearly, given (4.34), logarithmic singularities of the curvature invariant (4.33) can 
not be canceled within the supergravity approximation. 

5 Conclusion 

In this paper, following [25, 27, 28] we attempted to construct a string theory dual to 
strongly coupled conformal expanding plasmas in Bjorken regime [38]. In order to have 
computational control we truncated the full string theory to supergravity approxima- 
tion, and focused on the well-established examples of the gauge/string dualities: we 
considered M = 4 SYM [1] and super conformal Klebanov-Witten [51] gauge theories. 
We used non-singularity of the dual gravitational backgrounds as a guiding principle to 
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identify gauge theory operators that would develop a vacuum expectation value during 
boost-invariant expansion of the plasma. Truncation to a supergravity sector of the 
string theory (along with parity invariance in the Bjorken frame) severely restricts a 
set of such operators. In the case of the Af = 4 SYM, there are only two such gauge 
invariant operators, while for the Klebanov-Witten plasma one has an additional op- 
erator. We constructed supergravity dual as a late-time asymptotic expansion and 
demonstrated that the gravitational boundary stress energy tensor expectation value 
has exactly the same asymptotic late-time expansion as predicted by Miiller-Israel- 
Stewart theory of transient relativistic kinetic theory [47, 48] for the boost invariant 
expansion. As an impressive success of this approach, one recovers by requiring non- 
singularity of the background geometry at leading and first three subleading orders 8 
[25, 27, 28] the equation of state for the plasma, its shear viscosity and its relaxation 
time, in agreement with values extracted from the equilibrium higher point correlation 
functions. Unfortunately, we showed that logarithmic singularity in the background 
geometry can not be canceled within the supergravity approximation. Moreover, given 
that the singularities appear to persist in arbitrary high metric curvature invariants, we 
suspect that relaxing the constraint of parity invariance in the Bjorken regime would 
not help. Indeed, relaxing parity invariance would allow for only finite number of addi- 
tional (massive) supergravity modes, which, as an example of Klebanov-Witten plasma 
demonstrates, would only allow to cancel logarithmic singularities in finite number of 
additional metric curvature invariants. 

We would like to conclude with several speculations. 
■ It is possible that though the full asymptotic late-time expansion of strongly coupled 
expanding boost invariant plasma is ill-defined within supergravity approximation, the 
first couple orders can nonetheless we used to extract transport coefficients, and thus 
can be of value to RHIC (and future LHC) experiments. To this end we note a cu- 
riosity that the relaxation time of the M = 4 and the KW plasma came up to be the 
same. Moreover turning on the scalar fields did not modify the extracted relaxation 
time. This might point to the universality of the relaxation time, not too dissimilar 
to the universality of the shear viscosity in gauge theory plasma at (infinitely) strong 
't Hooft coupling (see [53] for a conjecture on a bound on relaxation times). Second, 
the relaxation time computed is the relaxation time for the equilibration of the shear 
modes. If there is substantial bulk viscosity at RHIC, one would also need an estimate 

8 At the third subleading order, logarithmic singularities of the metric curvature invariants remain. 
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for the corresponding relaxation time — this is where non-conformal gauge/string du- 
alities might be useful. 

■ It is possible that the singularity observed in the supergravity approximation is gen- 
uine, and can not be cured by string theory corrections. A prototypical example of 
this is the Klebanov-Tseytlin solution [54]. There, one does not expect a resolution 
of the singularity by string corrections (while preserving the chiral symmetry and su- 
persymmetry) simply because the corresponding (dual) gauge theory phase does not 
exist. By the same token, the observed singularity in the string theory dual to strongly 
coupled boost invariant expanding plasma might indicate that such a flow for a con- 
formal plasma is physically impossible to realize. For example, a Bjorken flow of a 
conformal plasma might always become turbulent at late times. In fact, at weak cou- 
pling, rapidity breaking instabilities were found [55, 56]. Although this was far from 
a hydrodynamical regime, the leftover logarithmic singularity appearing at the third 
order might be a manifestation of such an effect. 
We hope to address these issues in future work. 
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